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I consider the shape of a deformed elastic shell. Using the fact that the lowest-energy, small 
deformations are along infinitesimal isometries of the shell's mid-surface, I describe a class of weakly- 
stretching deformations for thin shells based on the Nambu-Goldstone modes associated with those 
isometries. The main result is an effective theory to describe the diffuse deformations of thin shells 
that incorporate stretching and bending energies. The theory recovers previous results for the 
propagation of a "pinch" on a cylinder. A cone, on the other hand, has two length scales governing 
the persistence of a pinch: one governing the relaxation of the pinch that scales with thickness as a 
— 1/2 power, and one that scales with thickness above which deformations again become isometric. 
These lengths meet at a critical thickness below which low energy deformations again become nearly 
isometric. 



As long ago as 1849, the Rev. Jellett pointed out that 
thin shells have a class of deformations that are privileged 
due to the large separation in scales between bending and 
stretching energy |l| . These isometries - deformations of 
the shell's mid-surface that do not stretch it - were the 
basis upon which Rayleigh's theory of small vibrations 
in shells was built [l, and were used more recently 
in a study of the rigidity of elastic torii 0, Q . In more 
modern applications, isometries appear in the guise of 
developable surfaces, characterizing all the forms that 
can be obtained by bending, but not stretching, a flat, 
elastic sheet 0. 

Under some types of forcing or constraints, however, an 
elastic sheet must stretch. The prototypical case is that 
of a crumpled sheet of paper , and more specifically an 
intrinsically-flat elastic sheet pushed into a cup to form 
a d-cone [8[. On a curved shell, matters are more com- 
plex. On the one hand, there are singular structures rem- 
iniscent of d-cones which sharpen with decreasing shell 
thickness [l3|- On the other, there are deformations of 
shells, such as "pinching" a straw, which result in very 
diffuse stretching over a large area [ll|. Finally, coexis- 
tence between diffuse and sharp features have been ob- 



served in wrinkled sheets 12-14 



This paper describes a generalization of shell isome- 
tries, analogous to Nambu-Goldstone modes in an or- 
dered phase. I derive a "universal" effective theory for 
the small- amplitude, diffuse (as opposed to sharp) defor- 
mations of thin shells. My approach is similar to that 
of Langeand Newell for the post buckling behavior of 



a shell [15|, [16[ , but applies more generally. The theory 
agrees well with previous results on the persistence of a 
"pinch" in a cylinder [llj , and generalizes easily to more 
general shells. This is illustrated by considering the de- 
formations of a cone, which qualitatively departs from 
those of a cylinder. 

The mid-surface of a thin elastic shell is represented 
by a set of coordinates {xi,X2) and functions X(a;i,a;2) 
describing its shape in space. For small strains, I follow 
standard practice in writing the energy as a sum of two 



terms, E = Est + Et [17|. The first term, due to in-plane 
stretching, is 



Est = 



J (fx VIA'J'^' 
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is the reference metric, g its 



2^g^^g^^ in terms of 



where t is the thickness, 
determinant, and A^^^^ - 
the two Lame constants. Alternatively, A = iV(l ^ ^'^) 
and 2/.J = Y/[l + i'), where Y is the Young's modulus and 
V the Poisson ratio. The strain is the difference between 
realized and reference metrics, — gij — cjij = diX. ■ 
i9jX — (jij. The second term in the energy, due to the 
bending of the sheet, is 



t 



hij][hki 



(2) 



where hij = N • didjX is the second fundamental form, 
N is the unit normal of the mid-surface, and hij is a 
reference bending form. When the reference metric can 
be realized as a surface immersed in three dimensions, a 
free elastic sheet will expel in-plane strain as t — >■ [l8( . 

The softest deformations of a shell are the isometries 
- deformations that change the shape of the mid-surface 
while leaving the metric invariant [l|. More precisely, I 
consider a shell with an iV-parameter family of isometries 
parameterized by e/, where / ranges from 1 to A^. An 




FIG. 1: (left) The coordinate r is the arc length on the un- 
deformed cone and 8 is the azimuthal angle, (middle) A cone 
pinched at the inside edge along an infinitesimal isometry. 
(right) A cone pinched at the inside edge along an isometry 
that decays along a Nambu-Goldstone mode. 
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order isometry requires d^\f=o{di'K ■ djX.) = for are 



any k < n. The larger n, the softer the deformation. In 
many cases, existence of a first- or second-order isometry 
is sufficient to imply an infinite order isometry [lij . 

Suppose we now consider the generalized class of de- 
formations in which itself can vary with position on 
the shell's mid-surface. Computing the in-plane stresses 
induced by such a deformation to first order, one obtains 



N 
1=1 



(3) 



where {T,y = d,iX ■ 9,X evaluated at = 0. Thus, 
if we restrict ourselves to this class of deformations, the 
stretching content of the deformation can be made van- 
ishingly small by making the gradient small. To this 
order, the total strain energy can be condensed into a 
simple form. 



; d 



\ V^gY^S^d^e'd.e-', (4) 



which must then be integrated over the area of the shell 



mid-surface. Notice, however, that the tensor SV, can 



generically depend on position. 

It is instructive to think of the example of an acous- 
tic phonon in a crystal which, like other Nambu- 
Goldstonc modes, is a gapless oscillation arising from 
a spontaneously-broken global symmetry. The broken 
translation symmetry of a crystal ensures that a long- 
wavelength oscillation of the atoms can have an arbitrar- 
ily low energy with increasing wavelength. Similarly, the 
ability to deform a shell along isometrics without stretch- 
ing energy cost ensures that the deformations resulting in 
Eq. ([3]) can also be made arbitrarily small. On the other 
hand, the stretching energy densities of non-isometric de- 
formations scale with Ytj'^. Thus, and in analogy with 
case of Nambu-Goldstone modes, I conjecture that gap- 
less, local isometrics dominate the diffuse, stretching de- 
formations shells. 

Unlike in the stretching energy, the bending energy 
will generically have terms depending on e^, die^ , and 
didje^ to quadratic order. However, because this pa- 
per is primarily concerned with very diffuse deformations 
and these terms are already down by compared to the 
stretching energy, I employ the approximation 



IJ 



(5) 



where Bjj is a position-dependent bending matrix. 
These two expressions for Est and Eb are the main result 
of this work. 

One can also compute the higher-order corrections to 
for the Nambu-Goldston-likc local isometrics. These 



„(2) 



e^d,e^d,iX ■ d^d^jX + e'dje'^d.jX ■ d,d,,X (6) 
+d^e' dje-' d,iX ■ d,.,X + e' e'^ d^,d^.j{d^X ■ djX). 



The last term of Eq. © vanishes for second order isome- 
trics but not necessarily in general. The first two terms 
can be compared directly to Eq. ([3]) and suggest that the 
limit of validity for the quadratic approximation used in 
Eq. dH) is, e-^d^jX ■ ^^^,,X < . Similarly, the third 
term of Eq. © yields the bound d^e-' d^iX ■ d,.iX < T/. 

What remains is to find expressions for the first-order 
isometrics. I require only that = 9g-r|o5ij ~ diX^ ■ 
djd.i |oXr+ a^Xo • \qX^. Therefore, D^T^ + D,Tj - 
2h^hij, where Di is the covariant derivative on the sur- 
face and /i-^ = N • d^j \ oX is the component of the isom- 
etry normal to the surface. Since only appears alge- 
braically, it can be eliminated to produce two equations 
in terms of the two unknown functions T/. 

I will focus on axisymmctric, closed shells, 
described by coordinates {^,0), X{r,9) = 

(p(r) COS0, pir) sine, d^^l - p'^o) ■ In these 
coordinates, r is the arc length along the symmetry axis 
and 6 the azimuthal angle around the symmetry axis. 
The isometrics can now be found in terms of a Fourier 
decomposition: (T;)^ = EmC^O^^e™^- Assuming that 
the can only vary in r, I obtain the equations 



Tr 



TV 



p (1 - p'2)3/2 
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(7) 



and 



pV^' 



(8) 



where the prime denotes a derivative with respect to r 
and I have suppressed extraneous indices. 

When m = 0, Eqs. [7] can be integrated in closed form: 
Tr = Vl-p'Hr)Ci, Tg = p^{r)C2, and h = p'{r)Ci, 
where Ci^2 are the two constants of integration. More 
generally, there is typically a two-parameter family of 
first-order isometrics for each Fourier mode (though this 
need not always be the case 0, S])- It is useful to asso- 
ciate these constants in terms of more physical boundary 
conditions. Denoting r = as the circle at which a defor- 
mation will be applied, I choose Tr{r = 0) = p{0)e^'"^^'^ 
and Tg{r = 0) = for one family of isometrics, and 
Tr^O and Tg{r = 0) = p{0)'^e'^'"'=° for the other. Thus, 



d.i.oX = 

9,2,0 X ~ 



1 



v^-p'{or 

pHr) 



Vl-p'(r)2a,X + p'(r)N 



p2(0) 



(9) 
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The first isometry controls translation along the sym- 
metry axis of the shell. The second isometry controls 
rotations about the symmetry axis. Similarly, the m = 1 
modes describe translations transverse to and rotations 
of the symmetry axis of the shell. The m = 2 mode on a 
cone is shown in Fig. [TJ 

The circle at r = plays a prominent role in Eq. (|9]) 
because of the choice of boundary conditions. For ex- 
ample, the m ~ 1 isometrics that rotate the shell do 
so about an axis perpendicular to the shell's symmetry 
axis and passing through the circle at r = 0. As a con- 
sequence, the expressions for Sfj and Bjj will depend 
explicitly on r, and would be superficially different were 
the boundary conditions applied elsewhere on the shell. 
However, the full symmetries of the shell energy under 
arbitrary rotations can be recovered by constructing lin- 
ear combinations of isometrics with different m. For ex- 
ample, combining a translation along the symmetry axis 
(m = 0) followed by a rotation of the symmetry axis 
about the origin (to = 1), one obtains an arbitrary in- 
finitesimal rotation. Similarly, translations and rotations 
can be combined with higher order modes (to > 2) to 
obtain deformations of the shell anywhere along the sym- 
metry axis. The choice of where to apply the boundary 
conditions is purely a convenience, because that is where 
we are imagining applying the deformation. 

Though there is no general analytic solution to Eq. 
([7]), the same analysis can be performed for each Fourier 
mode. For notational simplicity, we separate the sum 
over TO and let the index / be either 1 or 2. Then the 
strain components become 



ml 

Ire = ^6^™'(r)r/'"e 



(10) 



i.jnO 



and Jgg^O{e^). 

To sec how these pieces fit together, I will first recover 



the case of a cylinder deformed at r = [11| . for which 
p{r) = R is constant. There are two solutions for each 
TO, which we normalize with either = i?e^™, Tg ~ or 
T,. = 0,Tg = R^e^™. For the stretching energy, we need 
only 

' m'^r^/R^ + -imr/R: 

imr/R 1 
in Eq. (j4]). The bending energy is given by Eq. ([5]) with 

{X + 2^i){m^ - 1)2 



S7j = 



(11) 



Bij = 



i?2 

t^^+TOV/i?' 

im^r/ R 



2 im^r/R 
1 



(12) 




FIG. 2: Eigenvalues of the dimensionless stretching matrix for 
a "pinch" at r = on a cone [p{r) = R + a r] with a = 10""'^ 
for n = — 2.75. Arrow points toward increasing a. The 
result for a cylinder is shown as a thick line. 



The Eulcr-Lagrangc equations for the (e-*^"*, £2"*) are 



(SYjdre 



12 



= 0, 



(13) 



where repeated indices are summed. After some algebra, 
£2™ can be eliminated for £^™, leaving the single fourth 



order equation for e 



2 Am 



t^m^[m'^ - 1)2 
12i?6 



0, (14) 



where 0{t^) indicates a term of order that I will sub- 
sequently neglect. Setting to = 2, Eq. (fT4)) predicts 
a characteristic decay length A = (2\/3)^^/2^-i/2^3/2^ 
Denoting a persistence length £p = 27r/A, as in Ref. 



Ill ] , one obtains the correct scaling, the precession of the 
elliptical cross-section, and a prefactor ipt^^^R~^^^ = 
27r/v2%/3 ~ 3.38 that is not too far from the simula- 
tion value of, approximately, 4 ll|. Another difference 



is that the independence of the prefactor on Poisson ra- 
tio. However, even in the small deformation analysis of 
Mahadevan et al. [ll|, the dependence on Poisson ratio 



These expressions are only valid for to > 0. 



is weak. 

Eq. (jl3p can be easily generalized to other axisymmet- 
ric shapes and solved numerically. The simplifications 
that occur for cylinders leading to Eq. ([Ti)) do not oc- 
cur more generally, however. As a prototype of the more 
generic case, I consider a cone, for which p = R + a r, 
deformed at one edge located at r = 0. When a = 0, this 
shape recovers the cylinder previously discussed; when 
a ~ \, however, the surface is a flat sheet with a circu- 
lar hole of radius R at its center. We first consider the 
two eigenvalues of SYj/ipR^) for to = 2 and various a, 
plotted in Fig. [2l and the corresponding eigenvalues of 
Bi,jR^/{X + 2fi), plotted in Fig. [S] As a decreases, the 
eigenvalues approach those of the cylinder, a = 0, shown 
as thick lines. 

Since the largest eigenvalue of the stretching matrix 
grows with increasing distance from the pinch, r, I will 
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FIG. 3: Eigenvalues of the dimensionless bending matrix for 
a "pinch" at r = on a cone [p{r) — R-'r a r\ with a — 10""'^ 
for n = — 2.75. Arrows point toward increasing a. The 
result for a cylinder is shown as a thick line. 

logio(4^'^'/-R) i ncreasing « 
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FIG. 4: Persistence length of a "pinch" at r = i? on a cone 
[p{r) = R+a r] with a = lO""'' for n = 0-2.75. The cylinder 
is shown as a thick line. The critical point at tc is denoted by 
a circle. The dashed lines have slopes —1/2 (bottom) and +1 
(top) respectively. 

treat it as a constraint. Setting drt^ to vanish along the 
largest eigenvector of S/j, one obtains the approximation 
mdr€^ = aidrC^ for r ^ R. Since the smallest bending 
eigenvalue is already down by a factor oit^/12, it can be 
safely ignored. Consequently, only the smallest stretch- 
ing and largest bending eigenmodes compete. For small 
r, the stretching eigenvalue dominates and the shell de- 
forms along an isometry locally; at larger r, however, 
the largest bending eigenvalue may come to dominate 
and relax the deformation. Thus, I define the persistence 
length as the length at which the ratio of the smallest 
stretching and largest bending eigenvalues is one. For a 
general cone, there are two r at which this occurs, which 
I denote < £p\ both plotted in Fig. Uas a function 
of thickness for different values of a. It is only in the sin- 
gular limit a ^ 0, that a single length scale ^p"* ^ t^^^^ 
emerges. 

For r < £p\ the stretching energy dominates while 



for e);' < r < tp\ the bending energy dominates. For 

r > Vp \ the stretching energy dominates again. The 
multiple scales occurring in the cone have a natural in- 
terpretation: a deformation initially decays with a char- 
acteristic length scale but this decay halts beyond 
tp' as the stretching energy comes to dominate again. 
The separation between these length scales vanishes at a 
critical thickness, tc, below which the stretching energy 
always dominates. In this regime, diffuse deformations 
are essentially isometric. The critical thickness is not ap- 
parent on cylinders, since <c — ^ as a — >■ 0. On the other 
hand, when a approaches 1, so that the cone becomes 
nearly flat, the normal component of the isometrics di- 
verges with /i^" - 0[(1 - q;2)-i/2] It is clear that both 
bounds of validity emerging from Eq. ([6|) will be violated 
and nonlinear terms in the strain, which couple modes 
with different m's, will become relevant. 

Finally, one can ask when one should expect a descrip- 
tion of shell deformations in terms of local isometrics 
to fail. As suggested by a simple scaling argument, the 
stretching energy gap of non-isometric deformations may 
still be sufficiently large that the Nambu-Goldstonc-like 
local isometrics dominate even when the quadratic ap- 
proximation for the energy fails. That would suggest that 
the methods employed here are applicable quite generally 
to the deformations of shells. To know this for sure, how- 
ever, one would have to consider the size of these gaps, 
a more complex calculation not undertaken here. It is 
tempting to conjecture that a shell will crumple or form 
a sharp feature precisely as these non-isometric modes 
become important. This will be the subject of a future 
investigation. 
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